arXiv:1507.02652v3 [cond-mat.quant-gas] 6 Jan 2016 


Self-similar non-equilibrium dynamics of a many-body system 
with power-law interactions 

Ricardo Gutierrez, Juan P. Garrahan, and Igor Lesanovsky 
School of Physics and Astronomy, University of Nottingham, Nottingham, NG7 2RD, UK 

(Dated: January 7, 2016) 

The influence of power-law interactions on the dynamics of many-body systems far from equi¬ 
librium is much less explored than their effect on static and thermodynamic properties. To gain 
insight into this problem we introduce and analyze here an out-of-equilibrium deposition process in 
which the deposition rate of a given particle depends as a power-law on the distance to previously 
deposited particles. This model draws its relevance from recent experimental progress in the do¬ 
main of cold atomic gases which are studied in a setting where atoms that are excited to high-lying 
Rydberg states interact through power-law potentials that translate into power-law excitation rates. 
The out-of-equilibrium dynamics of this system turns out to be surprisingly rich. It features a self¬ 
similar evolution which leads to a characteristic power-law time dependence of observables such as 
the particle concentration, and results in a scale invariance of the structure factor. Our findings show 
that in dissipative Rydberg gases out of equilibrium the characteristic distance among excitations 
— often referred to as the blockade radius — is not a static but rather a dynamic quantity. 

PACS numbers: 67.85.-d, 05.30.-d, 32.80.Ee, 05.65.+b 


I. INTRODUCTION 

Self-similar behavior is ubiquitous in physics. It 
emerges for example in condensed matter systems near 
criticality [T], in biological systems [2] as well as in com¬ 
plex networks [3]. Systems with this property “look the 
same” under a rescaling of characteristic length and/or 
time scales. This property often entails a drastic re¬ 
duction of the details that are necessary for a theoret¬ 
ical description, such as specific parameter values or ini¬ 
tial conditions. Self-similar behavior can also occur in 
cold atomic ensembles, currently used as an experimental 
platform for exploring fundamental questions of far-from- 
equilibrium physics such as the approach of interacting 
many-body systems to (thermal) equilibrium [4]-[6]. For 
example, self-similar behavior is predicted for the relax¬ 
ation of a gas of quenched Bose-condensed atoms as a 
consequence of the existence of non-thermal fixed points 
EIH]. Moreover, it is argued [9] that many of these char¬ 
acteristics were also present during the evolution of the 
early universe leading to the idea of using cold atomic 
systems as analogue systems for addressing problems of 
relevance to cosmology or particle physics. 

In this work we introduce a simple far-from- 
equilibrium scenario in which a non-trivial relaxation dy¬ 
namics is driven by power-law interactions. The moti¬ 
vation stems from recent experiments in the domain of 
cold atom physics which explore the laser excitation of 
atoms to high-lying electronic states in which they in¬ 
teract with power-law potentials. While complementary, 
the scenario we study below has analogies with that of 
the above-mentioned work on Bose-condensed gases 13 - 
[9]. In particular, evolution which is initially fast and un¬ 
correlated is succeeded by slow and strongly correlated 
self-similar growth. We show that the physics of this 
dynamics is governed by an effective particle deposition 
process where the rates depend on the distance to other 


deposited particles as a power-law. The scale-invariant 
nature of the deposition dynamics is revealed analyti¬ 
cally and confirmed numerically by means of extensive 
numerical simulations. Moreover, we show that the re¬ 
laxation dynamics crosses over into a mean-field regime 
when the power-law exponent becomes equal to the sys¬ 
tem dimension. Beyond providing an understanding of 
non-equilibrium processes governed by power-law inter¬ 
actions our results reveal new insights into the relaxation 
behavior of gases of interacting Rydberg atoms, which are 
currently widely employed for the experimental study of 
many-body phenomena pTUlG] . Our results indicate that 
the characteristic minimal distance between Rydberg ex¬ 
citations — which is often referred to as blockade radius 
— is not generally a static quantity but can in the pres¬ 
ence of dissipation acquire a non-trivial scale-invariant 
t ime- dep endence. 


II. EFFECTIVE DYNAMICS 

The dynamical problem we study below emerges from 
considering a collection of N atoms on a cubic lattice in 
d dimensions. Each atom is described with two levels, a 
ground state and the excited Rydberg state, denoted as 
1^) and It) respectively. Atoms in the Rydberg state at 
positions Vk and interact through a power-law poten¬ 
tial = Ca/|r/c — with exponent a, while the 

much weaker interactions involving atoms in the ground 
state are neglected. Typically encountered exponents are 
Of = 6 (van der Waals interaction) and a = 3 (dipolar 
interaction) [17]. The atomic states |t) and ||) are res¬ 
onantly coupled by a laser field and at the same time 
subject to strong dissipation which leads to the rapid de¬ 
phasing of superposition states, e.g. |t) + |i). The evo¬ 
lution of the system is governed by a (classical) spin-flip 
dynamics in which a change in the state of site k, whether 
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FIG. 1. (Color online) Deposition process vs. full Rydberg gas dynamics, (a) Concentration c{t) of excitations as 
a function of time in a ID chain oi N — 10^ atoms interacting through a power-law potential {a — R — 15) for the full 
process (deposition and removal; blue dashed line) and for pure deposition (red continuous line). The black dashed segment 
corresponds to a power-law with exponent dl(2a + d) = 1/7. See text for details, (b) Same results for a iV = 100 x 100 2D 
square lattice. The black dashed segment corresponds to a power-law of exponent 1/4. (c) Natural logaritm of the deposition 
rate T^ [Eq. 0] as a function of ffc in a ID lattice {R — 15) with excitations being present on sites 1 and 20. (The most 
relevant cases, a = 3 and 6, are accentuated, and the order in which the lines appear starting from the top, a — and moving 
down to the bottom line, a = 6, is that of increasing a values.) (d) Sketch of the deposition model in a 2D lattice. Red (blue) 
discs represent excited (ground state) atoms. The darker the blue the higher the deposition rate for subsequent excitations. 
Encircled discs indicate the most likely positions for deposition. 


an excitation |^) ^ |t) or a de-excitation |t) ^ |i), oc¬ 
curs with an (operator-valued) rate (see Ref. [18]) 


(where jt) ^ |i) is not allowed) governed by the following 
master equation 


q-' = i 


n 2 


m^k 


T/c 


(1) 


Here R [T9| parameterizes the interaction strength, Uk = 
is the excitation number operator and Vk = r/^/a 
are position vectors in units of the lattice constant a. 
The term “ t*m|^ strongly correlates 

the atoms, giving rise to a slowing down of the dynamics 
in the vicinity (within a distance of the order of R) of an 
excited atom. In this sense, the dynamics is affected by 
kinetic constraints of the type that are usually considered 
in simplified models of the glass transition [20| . This de¬ 
scription of the Rydberg gas dynamics has been recently 
used in a number of theoretical works and was 

also successfully employed to model experiments nansj. 
A similar perturbative approach to the derivation of ef¬ 
fective dynamics in the limit of strong dissipation was 
proposed in Ref. m- 

We consider an ‘empty’ initial state, \C) = * * * i)- 

One key observation that greatly helps to simplify the 
following analysis is that excitation processes ||) ^ |t) 
strongly dominate over de-excitations |t) ^ |i), as will 
be shown below. Thus, the effective reversible dynamics 
is in essence a pure (and irreversible) deposition process 


dt\P{t)) = [4 - (1 -nfe)] \P{t)). (2) 

k 

The state is given by the vector \P{t)) = Xlc 
with P(C; t) denoting the probability of finding a specific 
configuration of excitations |C) at time t. The operator 
creates an excitation at site k, or — in the deposition 
picture — deposits a particle in site k. Figure shows 
the concentration of excitations c{t) = '^j^{'nk{t))/N as 
a function of time for the case of a = 3. We see from 
Figs. 0a) and[^b) [ID chain and a 2D square lattice, re¬ 
spectively] that irreversible deposition processes indeed 
dominate the growth regime, similar to what was sug¬ 
gested recently in Ref. [28| . Deviations become apparent 
in the long-time limit when the system approaches its sta¬ 
tionary state: for the reversible process this is the fully 
random state, limt^oo c(t) = 1/2, while for the deposi¬ 
tion process it is the (absorbing) state of a full lattice, 
limt^oo c(t) = 1. Note, that in the case of Rydberg gases 
this long-time limit will typically not be achieved due to 
the finite lifetime of excited atoms. However, while a fi¬ 
nite decay rate eventually makes the system settle into 
a nonequilibrium stationary state, the initial phases of 
the growth dynamics are well within the reach of current 
experiments, and they turn out to be virtually indistin- 
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guishable from those observed without decay, according 
to both numerical [18] and experimental evidence [29] , 


III. DEPOSITION MODEL 

In the strongly correlated regime, i.e. when the inter¬ 
particle distance becomes smaller than i?, the deposition 
process is amenable to an approximate analytical treat¬ 
ment. This is most transparently explained in ID: Con¬ 
sider two excitations in a ID chain separated by distance 

1. Then the next excitation will be deposited with very 
high likelihood at (essentially) a distance 1/2 from the 
two existing excitations, effectively rescaling the inter¬ 
excitation distance by the factor 6 = 2. The reason for 
this is that the deposition rate is highly peaked at the 
mid point between two excitations as shown in Fig. [^c). 
The process continues until the natural cutoff scale — 
the lattice spacing a — is reached, which makes further 
subdivision impossible. A similar process is expected to 
occur in 2D, where the same reasoning, starting for ex¬ 
ample from a square of excitations, can be applied with 
a rescaling constant b = ^/2 [see Fig. Bd)] . That the 
actual dynamics indeed follows a rescaling behavior is 
nicely observed in Fig. [^a) (details further below). 

For a more quantitative analysis we consider the distri¬ 
bution 7r(/,t) of distances I between nearest excitations. 
Formalizing the above considerations suggests that the 
dynamics of the deposition process is approximately de¬ 
scribed by the following master equation 

(hJ\d id 

dtAl, t) = i^r(M)7r(W, t) - t), (3) 

where 7r(/, t) is the probability distribution of the distance 
between nearest excitations I at time t. Here, d = 1, 

2, or 3 gives the dimensions of space and b = 2/^/d is 
the scaling parameter. The rate at which the rescaling 
step I ^ l/b occurs is r(/) « z~‘^R~‘^^P^(l){b) as soon 
as there are excitations within distances shorter than i?, 
according to Eq. (1). Here, 2 : is the coordination number 
of the lattice (z = 2d as we consider d-dimensional cubic 
lattices) and 0(6) is a geometric factor that depends on 
the space dimensionality. For instance, in ID 0(6) = 

where different terms in the 
sum correspond to different “shells” of excitations. We 
can therefore rewrite Eq. (§ as 

drn{l,T) = /j-j Tr{bl,T)-/j-j {i) 

where L = is the linear size of the lattice. In or- 

der to simplify the equation, we have absorbed the fac¬ 
tor z~‘^(j){b)L‘^'^^ which is fixed for a given lattice 
and interaction potential, into a rescaled time r. The 
fact that 0(6) does not play any role in the derivation 
of the self-similar behavior we report means that the ini¬ 
tial condition shown in Eig. (d) is simply a conve¬ 
niently simple example for the purpose of illustration. 


As the normalized inter-excitation distances l/L follow 
1 ^ 1/6^ 1/6^ ^ •••, we can rewrite Eq. Q in terms 
of 7r^(T), the probability of having an inter-excitation 
distance I = b~^L for n = 0,1, 2,..., 

drMr) = ( 5 ) 

where the rhs can be most economically expressed in 
terms of fn{r) = as fn-i{r) - fn{r). 


What we have in Equation ^ is the master equation 
in terms of 7r(/,t) when the deposition is such that I 
takes on a discrete set of values. Nevertheless, in gen¬ 
eral the initial condition 7r(/,0) is not a delta function, 
but a random configuration that arises from the initial 
creation of independent, distant excitations. Moreover, 
the fact that sometimes excitations do not occur exactly 
at the maximum of the rates forces us to consider con¬ 
tinuous distributions where the rescaling process with 
parameter 6 takes place simultaneously at slightly dif¬ 
ferent scales. To address these issues, we move on to 
the continuum, I = b~^L for x > 0, and Taylor ex¬ 
pand f{xR) (assuming analyticity with respect to x): 

f{x-l,T)-f{x,T) = -J2^^i^, In 

this setting, we can rewrite the master equation back in 
terms of I using the function n{l^T) = (//I/)*^^^+^^7r(/,r) 


^r7r(/,r) 


~ (_l)P+i 

^ p! 

p=l 


[(log (6) 


(6) 

Plugging in the wave-like ansatz 7r(/,t) = F{r — g{l)) 
(we expect the distribution to shift towards shorter and 
shorter distances as time goes by), we obtain 


(Z/L)-(2«+^) 


- log b - g{bl))blg^^\bl) (7) 

J^F(^\T-g{bl))blg^^\bl) 

- g{bl)) hlg^'^\bl) 
i"(i)(T - g{bl)){blfg^^\bl) + • • • 


stopping at order p = 2. Eor higher orders on the rhs 
one obtains more terms that are (up to some constant 
factor, including functions of 6) products of F(-) or its 
derivatives and {bl)^g^^\bl) raised to a certain 

power, where g^^\bl) is the m-th derivative of g{') eval¬ 
uated at bl. As the I dependence has to cancel on both 
sides for the equation to hold for all times, we are left 
with g{l) = and therefore (after absorbing L 

into the undetermined function F{')) 


n{l,T) = r(2«+‘i)F(r - r(2«+‘')). 


(8) 


The mean inter-excitation distance as a function of 
time, (/(r)) = j (i//7r(/,r)/ j (i/7r(/,r), is, after chang¬ 
ing the integration variables to ^ = r — /-(2a+c^)^ 


{l{r)) 


_ 1 

— p- 2o!.-\-d 


j dy{l- y/r) 
f dy(l- yjx) 


2a + d F(^) 

F[\j) 


(9) 
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FIG. 2. (Color online) Scale invariant evolution, (a) Ex¬ 
citations (blue dots) in a 2D square lattice of = 256 x 256 
atoms, 0 = 3 and R = 15. The first column shows a 
subsystem of size h x h sites for /i = 128 at three dif¬ 
ferent times: ti = 1, t 2 = 10^, ts = 10®. The second 
column shows the configuration within a box of linear size 
h = — 54 at times t 2 and ts. The third 

column shows the configuration within a box of linear size 
h = — 23 at time ts- (b) Snapshots of the 

structure factor S{k) taken at different times, (c) S{k) as 
a function of the scaled wavenumber k/^/c(T). All the lines 
shown in panel (b) collapse onto a single curve. 


For any dilute initial condition, F(-) is non-zero for small 
values of y. On the other hand, the process spans sev¬ 
eral orders of magnitude in time (see Fig. 0 (a, b)), so 
for sufficiently long times there has to be a strong can¬ 
cellation in the argument of F, and must grow 

essentially in parallel with r. This means that eventually 
(1 — yI t) ~ 1 and, writing back in terms of the original 
time units, (/(t)) ^ or, equivalently, 

c{t) = ^ (10) 

This result is in agreement with the prediction obtained 
and verified in Ref. [18] by a more restrictive reason¬ 
ing that did not consider the time dependence of the 
full distribution 7r(/,t). Eq. ( [1Q| ) justifies the use of the 
exponents 1/7 and 1/4 in Fig. (a, b). By the same 
reasoning we can rewrite 7r(/,t) as given in Eq. (§ as a 
scale invariant function for long times: 

Tr{i, t )« r(11) 

showing that the full distribution of distances to nearest 
excitations evolves in a self-similar way. Physically, this 


means that the system looks the same at any time, only 
the characteristic distances become shorter and shorter. 
The prediction of self-similarity from the approximate de¬ 
scription in terms of 7r(/,t) appears robust for the actual 
deposition problem, as seen from kinetic Monte-Carlo 
simulations of Eq. (§: Fig. I^a) shows that the con¬ 
figurations of a 2D system at different times, and thus 
very different concentration, look statistically similar if 
lengths are scaled according to Eq. The data in 

Figs. (b,c) confirm this observation through the col¬ 
lapse of the structure factors at different times. 



FIG. 3. (Golor online) Scale invariance of 7r(/, t) in ID and 
2D for 0 = 3. (a) Distribution of distances between nearest 
excitations 7r(/, t) and rescaling according to Eq. ( |11| for a — 
3 in a ID chain (iV = 2^®, R = 15). Lines of different color 
correspond to different times. The ins et s hows the collapse 
of the distribution functions under Eq. 0 [distributions are 
shown as a function of I , and multiplied by ^-1/(2“+^) 

to account for the change in the measure of the distribution]. 
(b) Same as panel (a) for 2D (N = 256 x 256, R = 15). 


Furthermore, the precise scaling relation in Eq.(ll) is 
also verified in the simulations. The numerically obtained 
distributions of distances between nearest excitations are 
shown in Fig. The insets show the rescaled distribu¬ 
tions, which nicely collapse onto a master curve. This 
works quite well even for small values of /, where the dis¬ 
creteness of the lattice introduces some roughness in the 
curves, especially in 2D [Fig. i(b)]. 


IV. DYNAMIC RYDBERG BLOCKADE 

The form of the rates Q leads to pronounced spatial 
anti-correlations, i.e. there exists a “correlation hole” 
between excitations, cf. Figs. Bd) and|^a). In the con¬ 
text of Rydberg gases the linear size of this correlation 
hole is often referred to as the blockade radius Rb • Since 
the concentration increases in time Rb will become a 
dynamic quantity and thus we encounter a dynamic Ry¬ 
dberg blockade. The time-dependence of Rb can actually 
be inferred from a macroscopic measurement, namely the 
analysis of the system-size dependence of the fluctuations 
in the number of excitations. The excitation number vari¬ 
ance at a given time in the non-equilibrium evolution of 
a system of size N is var(A^) = {n{N)‘^) — (n(A^))^, where 
n{N) = After some algebra, this can be written 
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in terms of 5ni = rii — c (as usual c = (n^)) as follows 


var(Ar) = N 





( 12 ) 


The anti-correlations due to the blockade effect are ex¬ 
pected to decay at a finite distance for sufficiently large 
N. Indeed, starting from a very small N (within the 
blockade radius), ^kj^Q{Sni6ni-^k) is expected to be¬ 
come more and more negative as N grows, until N 
reaches beyond the Rydberg blockade radius, Srii and 
for sufficiently large If^—become independent 
{{6ni6ni-^k) = 0 ), and therefore satu¬ 

rates. In Fig. [^we show numerically obtained Y8iT{N)/N 
as a function of A/" in a ID lattice. The fact fhat the 
curves become flat for sufficiently large N is in agree¬ 
ment with the simple reasoning enunciated above. As 
expected, for small sizes a more complex dependence of 
the fluctuations on N is observed. 


atom can be located at any distance from an excitation). 
When N > 2Rb, we obtain the aforementioned satura¬ 
tion '^k^QiSriiSrii-^k) = (2Rb ~ 1)^(1 — c)(—c). Summa¬ 
rizing, 

iVc(l-c){l-(iV-l)c}, iV<2i?B 
'' ^ \ A^c (1 - c) {1 - {2Rb - 1) c} , N>2Rb. 

(13) 

Despite the simplicity of the approximation, Eq. is 
shown to be essentially valid in Fig. where we show 
var(A/') for different times together with the approxima¬ 
tion ( p!^ . The data collapses upon rescaling the system 
size by Rb oc t~ highlighting its usefulness for study¬ 
ing the dynamical nature of the Rydberg blockade from 
macroscopic observations. The left hand side is calcu¬ 
lated numerically (or it is measured in an experiment), 
which allows us to extract the value of Rb- 


V. MEAN-FIELD DYNAMICS ARISING FROM 
LONG-RANGE INTERACTIONS 



FIG. 4. (Color online) Excitation number variance as 
a function of the system size. Variance in a chain of 
2^® atoms (a = 3, R = 15) for (sub-)system sizes N = 
2,4, ...,2^^ (averaged over 1000 realizations). The dashed 
curves correspond to the expression in Eq. ( |13[ ) applied to 
each c, where Rb is extracted from the largest N consid¬ 
ered. Different lines correspond to concentrations starting 
from approximately c = 0.09 (light green line at the bot¬ 
tom) to c = 0.30 (pink line on top) [lines have been displaced 
vertically to improve the visibility]. Inset: Number variance 
curves normalized byAc(l —c){l — (2Rb — 1) c} and plotted 
as functions of N/Rb- 


To understand the small N behavior, we consider a 
hard-objects description |3Q1433] . i.e. we assume that 
there are no two excitations within a radius Rb- For 
the sake of simplicity, we focus on the ID case, al¬ 
though the same reasoning is expected to be relevant 
in lattices of higher dimensions. If A" < 2Rb, then 
^kfoiSriidrii^k) = {N - l)c(l - c)(-c). Here, c is the 
prob^ability of having = 1 , i.e. an excitation at a 
generic site i, (1 — c)(—c) is the product of the fluctu¬ 
ations at i (It)) and i k (|t)), and N — 1 results from 
the summation. We do not consider the case rii = 0 be¬ 
cause then {6ni6ni-^k) = [c(l — c) + (1 — c)(—c)] = 0 , 
where the term in square brackets includes the fluctu¬ 
ation for It) and |t) states weighted by their probabil¬ 
ity of occurrence ai i k (considering a ground state 


Finally, we consider the dependence of the deposition 
dynamics on the power-law exponent a. Qualitatively 
there are two different situations, that of short-range in¬ 
teractions {a > d) and that of long-range interactions 
(a < d) |34]. Figure shows the time evolution of 
the concentration in both cases. Whenever the power- 
law rates are short-ranged, the growth is exponential 
for very short times, and then follows the scaling law 
c{t) = cf. Eq. ( [Tq| ). In contrast, in the long- 

ranged case the concentration grows in an a-independent 
manner, c{t) ^ H, with 7 1/3 throughout (or, more 

precisely, after an extremely short exponential growth 
whose duration decreases with growing system size). The 
behavior of the concentration for a < d can indeed be un¬ 
derstood by a mean-field analysis, as we presently aim to 
show. 

The deposition dynamics is governed by Eq. ([^, with 
rates given by Eq. 0 - The average of the number oper¬ 
ator corresponding to site j, {nj{t)) = {—\nj\P{t)), with 
|—) = X]^|C), evolves in time according to 


dt{nj{t))= E(-mni [o-+ - (1 - rik)] \P{t)) (14) 


= (-|r 


J 


7-(l-nj) \P{t)) 


= (-|B(l-n,-)|P(i)) 


where the mean-field lack of correlations between sites is 
exploited. Moreover, given the equivalence of all sites, 
nj{t) dtc{t) = (—|r(n(t)) (I — n(t))|P(t)), where 

as usual c{t) = {n{t)). We can approximate this equation 
by replacing expectation values of products of operators 
by products of expectation values of operators: dtc{t) ^ 
r(c(t))(I — c(t)). The resulting mean-field equation for 
the time evolution of the density is 


dtc{t) = 


1 - c{t) 


(15) 
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Here, the geometric factor \^k\~^ converges for 

a > d and grows unboundedly for a < d as is in¬ 
creased. It can be shown that the system asymptotically 
reaches the stationary solution = 1 . 

For Of > d, we expect the existence of an initial stage, 
c{t) ^ 1, for which the denominator in Eq. © is essen¬ 
tially one, and there is an initial exponential growth 

c{t) = 1 — exp(—t). (16) 


ation 

c{t)^l-exp{-t/{l + F^R^^)). (18) 

So the final stages of the mean field dynamics are ruled 
by exponential growth in all cases, as seen in Fig. 

To study the intermediate time regime, for which c(t) 
is still considerably smaller than 1 , it is useful to expand 
the implicit solution in powers of c{t): 


Physically, this corresponds to the creation of indepen¬ 
dent, distant excitations that eventually become the ini¬ 
tial seed for the deposition process. In the case of a < d, 
due to the unbounded growth of |r/c|~^ with 

the system size, this regime is expected to be negligibly 
short for sufficiently large N. This is in agreement with 
the results reported in Fig. 



FIG. 5. (Color online) Time evolution of the concentra¬ 
tion for different power-law exponents a. (a) Concen¬ 
tration c{t) in a square lattice {N = 317 x 317, R = 15) with 
long-range interactions, a < d — 2 (averages of 30 realiza¬ 
tions). The larger the value of o, the faster the dynamics 
(and therefore the higher the position at which the line ap¬ 
pears in the plot). The dashed segments superimposed to 
each line show the result of a power-law fit c(t) ~ F. (b) 
Same as in (a) for a > d. In this case, the larger the value 
of o, the slower the dynamics (and therefore the lower the 
position at which the line appears in the plot), (c) Growth 
exponent 7 vs. interaction exponent a in the same 2D square 
lattice (red squares) and a ID chain (black dots; N = 10^, 
R = 15). Dashed lines correspond to 7 = 1/3 (mean-field 
exponent, for a < d) and ^ — dj ( 2 a -h d) (valid for a > d). 


t = c(i)+c(t)V2+(l + F2i?2«) [c(t)3/3 + 0(c4)] . (19) 

If [FaR'^c{t)]‘^ ^ 1, we obtain the following algebraic 
growth: 


c{t) « [ 3 /(F 2 i 72 «)]V 3 ^i/ 3 ^ 


( 20 ) 


This is in agreement with the exponent 7 ~ 1/3 observed 
for a < d in Fig. as is the fact that the prefactor of the 
power law increases as a becomes larger. For long-range 
interactions, a < d, the dynamics is indeed effectively of 
the mean field type. 

In conclusion, there are three distinct dynamical 
regimes qualitatively identical to those reported in Ref. 
[18]. For very short times, the creation of independent, 
distant excitations makes the concentration grow expo¬ 
nentially in time (for large N this is only observable if 
a > d), as predicted at the mean field level. For very 
long times, the system reaches exponentially the station¬ 
ary state. Throughout most of the non-equilibrium evo¬ 
lution, the concentration grows algebraically according 
to the self-similar dynamics, Eq. (10), for a > d, and 
following a mean-field dynamics with exponent 7 = 1/3 
in the long-range interacting case. 

To gain more insight into the dynamics of the mean- 
field regime, we study the distribution of interparticle 
distances 7 r(/,t) for a < d. In the absence of spatial 
correlations, which is the characteristic of mean field ap¬ 
proaches, the probability of finding zero excitations in a 
small spherical volume 6V (in units given by a^) at time 
t is P{0;6V) = 1 - c{t)6V + 0{SV‘^). Therefore, the 
probability of finding no excitations in a volume V + SV 
is 


P(0; Vf6V) = P(0; V){1 - c{t)6V + 0{6V^)), ( 21 ) 

and thus (P(0; V F 6V) - P(0; V))/SV = -c{t). In the 
limit of 5V ^ 0 , we obtain the normalized distribution 

P(0; V) dV = dV. ( 22 ) 


For a better understanding of the other regimes in¬ 
volved, we obtain the following implicit solution to Eq. 
(15) by separation of variables: 


rp2 D2a 

( = -(1 + log (1 - c(t)) - (2c{() + c{tf) ^ 

(17) 

For long times, c{t) is approximately 1 and the logarith¬ 
mic term is dominant, which gives an exponential relax¬ 


So far in this paragraph, the reasoning has followed the 
standard Poisson process derivation. As F = A(// 2 )^ 
(where A is a geometric factor), we write in terms of /, 
and using the previously established notation, 

tt{ 1, t)dl = d{A/2'^)c{t)l'^-^e-<*^^^/^^^^^dl. (23) 


The time dependence can be made explicit by using the 
mean-field result displayed in Eq. (20), c{t) ~ ^ 
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where for convenience we define B = 
resulting in 

nil, t)dl = (24) 

In terms of a rescaled variable x = the distribution 

can be written as 

7r(x, dx, (25) 

where we have explicitly included the dimensions of space 
as a superscript in It turns out that the full 

distribution is also scale invariant in time in the long- 
range interacting case, and a rescaling similar to that seen 
in the case a > d^ but with exponent l/Sd, should thus 
be possible. This is indeed confirmed in Fig. which is 
analogous to Fig. but for a = 0.8 instead of 3. Due 
to the oscillations observed in ID at small distances, we 
limit ourselves to relatively early times in that case, to 
make the approximate adequacy of the collapse at least 
partly visible. These oscillations seem to indicate that 
the mean-field assumption does not hold perfectly, and 
some spatial correlations are still visible at a = 0.8 at 
the level of the full distribution function. 



FIG. 6. (Color online) Scale invariance of 7r(/, t) in ID and 
2D for a = 0.8. (a) Distribution of distances between nearest 
excitations 7r(/, t) and rescaling according to Eq. ( |20| ) for a = 
0.8 in a ID chain {N = 2^®, R = 15). Lines of different color 
correspond to different times. The inset show s the collapse 
of the distribution functions under Eq. ( [^ [distributions 
are shown as a function of 1 and multiplied by to 

account for the change in the measure of the distribution. The 
black line corresponds to the analytical prediction 7r(x, = 

d{A/2^)Bx^~^e~^^^‘^ (see text for an explanation)], (b) 

Same as panel (a) for 2D {N = 256 x 256, 15). 


Besides knowing that there has to be an approximate 
scale invariance for a < d due to the approximate validity 
of Eq. (25), we can use that equation to calculate explic¬ 
itly the functional form of the rescaled distribution. In a 
one dimensional chain, we obtain 7r(x,t)^^^ = 
while in a two-dimensional lattice the distribution is 
7r(x,t)^^^ = 27rBxe~^^^ . In these expressions, the sym¬ 
bol TT without explicit dependence on any variables is just 
the constant tt (the ratio of a circle’s circumference to its 
diameter), and it should not be confused with the distri¬ 
bution 7r(x, The corresponding curves are shown as 
continuous black lines in the insets of Fig. (a,b). We 
want to emphasize that no fitting procedure has been 
used, as the functional form is completely closed except 


for the parameter 5, which can be obtained numeri¬ 
cally for a lattice of a given topology and size. In fact, 
B changes depending on the dimensionality because 
does, and, for a cubic lattice oi N = 2^^ sites, its value is 
B = 0.0198 for d = 1 and B = 0.00214 for d = 2. While 
in the case of d = 1 [Fig. |^(a)] the agreement is not very 
good (which is to be expected, as it presupposes a con¬ 
stant growth as I becomes smaller that is not realistic in 
the presence of blockade effects -even in this largely, but 
not perfectly, mean field regime), in d = 2 [Fig. |6](b)] 
the fit is fairly good. So it is clear that the distribution 
of inter-excitation distances in the long-range interacting 
regime also shows at least approximate scale invariance, 
and that the distribution is not far from that one would 
find in the idealized mean-field case, i.e. in the complete 
absence of correlations. 


VI. CONCLUSIONS 

We have introduced a simple far-from-equilibrium sce¬ 
nario in which a non-trivial relaxation dynamics is driven 
by power-law interactions that result in rates that depend 
on the distance to existing excitations as a power-law. 
The motivation originates from the study of Rydberg 
gases far from equilibrium. We provide a simple model 
for the deposition process that captures the essential 
physics. The scale-invariance of the deposition dynamics 
is revealed analytically and confirmed by means of exten¬ 
sive numerical simulations of the original problem (with¬ 
out de-excitations). Our results indicate that the block¬ 
ade radius acquires a scale-invariant time-dependence in 
the presence of dissipation. Moreover, we study how the 
dynamics depends on the exponent of the interactions, 
and show a crossover into a mean-field regime when the 
power-law exponent becomes equal to the dimensions of 
the lattice or smaller. 

Despite the evident differences in the underlying physi¬ 
cal models, it is worthy of mention that a scale-invariant 
dynamics has also been theoretically established in ob¬ 
servables characterizing domain growth and phase or¬ 
dering in classical models of magnetism that have been 
quenched from a disordered phase to an ordered phase, 
even in the absence of power-law interactions [35]. In- 
triguingly, this is another setting in which the dynamics 
can be thought of as a deposition process. It is not clear 
to us whether deeper relations between that problem and 
the one that concerns us here exist (though the fact that 
the power-law exponents differ speaks against the use of 
conventional universality arguments). 

Whether connections to problems studied in the past 
are eventually found or not, the presented out-of- 
equilibrium setting is probably one of the simplest man¬ 
ifestations of a many-body evolution governed by power- 
law rates, and therefore we expect that a similar dy¬ 
namics may be observed in systems of quite a differ¬ 
ent nature than the dissipative Rydberg gases that in¬ 
spired its study. In view of the experimental realization 
















through Rydberg atoms it is interesting to ask how much 
of the observed features actually persists in the quantum 
regime, i.e. when the strong noise condition is lifted. 
While computationally unfeasible this could very well be 
explored in the most recent generation of experiments, 
e.g. pm [Ml [36] . 
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